Connecting wave functions at a three-leg junction of one-dimensional channels 
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We propose a scheme to connect the wave functions on different one-dimensional branches of 
a three- leg junction (Y-junction). Our scheme differs from that due to Griffith [Trans. Faraday 
Soc. 49, 345 (f953)] in the respect that ours can model the difference in the widths of the quasi- 
one-dimensional channels in different systems. We test our scheme by comparing results from a 
doubly-connected one-dimensional system and a related quasi-one-dimensional system, and we find 
a good agreement. Therefore our scheme may be useful in the construction of one-dimensional 
effective theories out of (multiply-connected) quasi-one-dimensional systems. 
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I. INTRODUCTION 

For a system which comprises quasi-one-dimensional 
(QID) channels, when only the low-energy regime at near 
the first subband bottom is considered, it can usually be 
modeled by a one-dimensional (ID) system. When the 
system is multiply-connected and consists of multi-leg 
junctions, the wave functions on the branches are usually 
connected at the junctions by the Griffith schemejiiSi^ 
the Shapiro scheme^^i^ or schemes alike. Since such for- 
mulations greatly reduce the calculational effort of com- 
plicated multiply-connected mesoscopic systems, they 
have been use d widely in the literat ure. F or example, 
see Refs. IllTollTlllTalTlTilTlll^ and the ref- 

erences therein. However, arguments which lead to these 
connecting schemes are kinematicalji*2ii^iSiSii and it is 
not clear what kind of junction in practice they de- 
scribe. Moreover, a comparison between the results of 
these schemes and that of the exact calculation of QID 
systems has never been done. It is the purpose of this 
paper to make a comparison between the Griffith result, 
the QID result, and the result due to a scheme we pro- 
pose in this paper. We find that for clean junctions of 
QID channels, the Griffith result is not even qualitatively 
in accord with the exact result. The scheme we derive 
gives a result that compares much better with the exact 
result. 

At a A^-leg junction of ID channels, the wave function 
continuity condition is a requirement that must be re- 
spected. Besides, the Griffith schemei*2i2i4 demands that 
the sum of the derivatives of the wave functions on the 
different branches at the junction is zero, i.e.. 



0, 



(1) 



where the directions of the coordinates are defined 
either diverging from or converging to the junction. 
This is the simplest way to impose the unitarity con- 
dition of no net current flows into the junction, i.e., 
^jRc 'il,'*{—idijji/dxi) = 0. When there is a magnetic 
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FIG. 1: (a) The original Y-junction of QID channels consid- 
ered in this paper, (b) the reformulated Y-junction of tight- 
binding channels, and (c) the effective Y-junction of ID chan- 
nels affiliated with the Griffith or our connection scheme. 



field, the requirement is rephrased as the sum of the co- 
variant derivatives is zero, i.e., ^i{d / dxi — ieJ^l)'ipi = 0, 

where a'' is the component of the vector potential par- 
allel to branch i at the junction. On the other hand, the 
Shapiro schemo^*^ directly demands that the scatter- 
ing matrix connecting the in-going and out-going waves 
at the junction be unitary, and a general matrix with 
free parameters is written down. When the spin degree 
of freedom is considered, these schemes are straightfor- 
wardly applied to each spin channel ii^*iSiiLi2iiS These 
schemes and the likes have been taken for granted and 
used widely in the literature. 



II. FORMULATIONS AND MODELS 

We approach the problem from another point of 
view. For a QID system with equal-width channels (the 
"width" is an ill-defined quantity in snaking channels but 
nevertheless we may talk about it when the curvatures 



2 



are small enough), we may approximate a three- leg junc- 
tion (Y-junction) and its branches [e.g., see Fig.^a)] by 
a tight-binding (TB) model as shown in Fig. ^b)i2£ 

The tight-binding model is described by a second quan- 
tised Hamiltonian 

H^J2clh,iCi, (2) 

where q is the annihilation operator of a spinless particle 
on site i, and /ly is a matrix element which is complex in 
general. The element /ly is called a hopping when i 7^ j, 
and an onsite potential when i = j. Defining a basis 
set {|i)} by |i) = c[|0), where |0) is the no-particle state, 
one can write the time-independent Schrodinger equation 
Hlip) = E\tp), where E is the energy of the particle, into 
the form 

J2ihij-E6,i)i;i^0, (3) 
j 

where 'ipj = is the TB wave function at site j. 

We define that the hopping exists only between nearest- 
neighbor sites, and be denoted by —t. The onsite poten- 
tial at site i is denoted by Vi + 2t. 

The magnitude of the hopping —t is obtained by the 
following argument. Let a QID channel be approximated 
by a finite-difference square grid, with three grid-points 
across the channel, one at the center and each edge. Then 
the distance between the grid-points will be m;/2, where w 
is the width of the channel, and the hopping in the finite- 
difference time-independent Schrodinger equation2i will 
be -t -2h^/{mw'^). We wiU assume the same hopping 
in our TB formulation. 

Away from the junction, the TB time-independent 
Schrodinger equation readsSi 

- ^(VJi+i - i^i) + ti-i,; - ^/^i.i) + {Vi - E)i>, = 0, (4) 

where E is the energy. In the long-wavelength 
limit it reduces, as it should, to the ID second or- 
der differential time-independent Schrodinger equation, 
-[h^/{2m)]dl^p{x) + [V{x) - E]-<p{x) = 0. 

At a Y-junction, the TB time-independent Schrodinger 
equation reads, 

(Vii - i^o) + W"2 - V^o) + (^3 - V^o) + ^ ~ ^° ^ ^ ^"0 = 0,(5) 

where the subscript "0" denotes the site at the junction, 
and "1" , "2" , and "3" denote the sites on the branches 
nearest to the site at the junction [i.e., in Eq. take 
i = 0, and j = 0, 1, 2, and 3]. It is seen that the Griffith 
scheme formulated in Eq. (QJ is recovered only when E — 
Vo + t = Q sX the junction. It is reasonable to set iJ = 
here since we are considering energies at near the band- 
bottom and E <^ t. But one still requires Vq = t to 
send the last term in Eq. Q to zero. In other words, 
the Griffith connection schemaiiSiii actually describes a 



Y-junction of QID channels with a repulsive potential 
with a strength of the order of t. Whereas in this paper 
we propose a connection scheme in the long-wavelength 
limit for a cZean Y-junction [i.e., Vq = in Eq. (Q] of QID 
channels. At a Y-junction of ID channels [see Fig. ^c)], 
we propose 

where the directions of the coordinates are defined to 
be diverging from the junction. If Eq. ^ is reached by 
dividing Eq. (jSJ by u'/2 and letting w — > (remember 
that E/t and Vo/t have been set to zero), the factor i' 
will be equal to 1. Adopting i> ~ I indeed results in a 
good enough qualitative comparison with the QID result. 
But we will see that choosing ~ 1.9 may bring the ID 
and QID results to a semi-quantitative agreement, which 
means that the term has been underestimated. The TB 
argument serves to bring out the 1/w dependence of the 
term, and the fixing of v will be discussed using concrete 
examples. The effect of the channel width is hence in- 
cluded, in contrast to the Griffith scheme (the case of 
= 0). The 1/w dependence results in an effect that is 
more prominent at smaller channel widths, and this un- 
derstanding may also help to relate studies on the quan- 
tum graph theorji^ to the practical experiments. The 
case of a general A^-leg junction can also be worked out 
likewise. 

In this paper we compare the Griffith and our schemes 
with the exact QID calculation in a chosen type of sys- 
tem. We calculate the transmission probability for a ID 
ring connected to two leads [see Fig. IS^a)], which is 
the simplest multiply-connected ID system, using the 
Griffith and our connection schemes at the Y-junctions. 
In addition, we also calculate the transmission proba- 
bility for a similar system, an annulus connected to two 
QID leads [see Fig.|2Ib)], using the exact mode-matching 
method. The two three-leg junctions in Fig. |2Ib) resem- 
ble the one in Fig.^a). Note that the transmission prob- 
ability is directly related to the experimentally measur- 
able conductanceiSi We will sketch how we have done the 
calculations, and we refer the readers to the literatures 
for more details. 

In a ID model as shown in Fig. El^a), the wave function 
on each line segment at a given positive energy E is a 
superposition of forward and backward traveling waves, 
i.e., 

i;, {x^) = A.e'^''^ + B,e-'''''\ i = 0, 1, 2, and 3, (7) 

where k = y2mE / h and m is the effective mass of the 
particle. The wavelength A is given by A = 27r/fc. The 
2^0.1,2.3 are the coordinates on the line segments corre- 
spondingly, and the coordinates have positive directions 
as that defined in Fig. |21 We define xi = X2 = X3 = 
at the right junction, and xq = 0, xi = Li, a;2 = L2 at 
the left junction, where Li_2 are the lengths of the arms 
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FIG. 2: (Color online) The doubly-connected ID system and 
the related QID system we consider, (a) The ID system is a 
ring (with arms labeled by 1 and 2) connected to two leads 
(labeled by and 3). The coodinate system Xi is defined 
for the line segment labeled by i (i = 0, 1, 2, and 3). The 
arrows denote the positive directions of the coordinates, the 
right Y-junction is defined at x\ = X2 = = 0, and the left 
Y-junction at xq = G, x\ = L\, and X2 =1/2. (b) The QID 
system is an annulus with two radially connected leads. Both 
the annulus and the leads have the same width. 



between the junctions [see Fig. Efa)]. The Ai (Bi) is the 
coefficient of a forward (backward) travehng wave. Since 
we consider particles incident from the left, we set = 1 
and i?3 = 0. Then the continuity requirement 



V'oUo=0 = -011:11=^1 = "021 



2\X2 = L2 



(8) 



and 



V'lUi=0 = l/'2U2=0 = -031:13=0, (9) 

and the Griffith unitarity imposition [following Eq. (^] 
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constitute an equation set which contains six equations 
with the six unknowns {Bq] Ai, Bi; A2, B2', A3} which 
have been defined in Eq. Q . Hence the transmission am- 
plitude A3 can be solved, and the transmission probabil- 
ity T = 1^3 [^ be found. We may also replace the Griffith 
unitarity condition by our unitarity condition [following 



Eq. ®] 

dxo 
and 
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and too the transmission probability can be solved. We 
will discuss the fixing of v later in this paper. 

Besides the mentioned ID model, we also solve a 
related QID model in a way as that of Xia and Li 
in Rcf. (22j. Consider an annulus with an inner and 
an outer radii of R ~ w/2 and R + w/2 respectively, 
and two leads of width if radially connected to it as 
shown in Fig. |2Ib). The wave function is governed by 
the two-dimensional (2D) differential time-independent 
Schrodinger equation. In a lead it can be expanded 
in terms of transverse modes (subbands) and longitu- 
dinal forward and backward modes, i.e., '!/'icad(2;, y) = 
^fLii'^i^^'^'^ + 6;e~*'^'^)sin(Z7ry/ii;), where x and y are 
respectively the longitudinal and transverse coordinates 
for the lead. The fc; and I are related by kf + (Itt/w)'^ = 
2mE/h^, where E is the energy (positive) of the par- 
ticle, and ki can be real or imaginary. In the annulus 
the wave function can be expanded by radial and angu- 
lar modes, i.e., Vannuius 6*) = Y^ii^M <l>i{kr)e'''^^ , where 
a radial mode is given by (t>i{kr) = ci Ji{kr) + diYi{kr), 
and k = \/2mE/h. The r and 6 are the radial and an- 
gular coordinates respectively; and the J; and Yi are 
the Bessel functions of the first and second kinds re- 
spectively. We demand 4'i\r=R-w/2 = for any 9, 
V'annuius|T-=ii-i-i«/2 = wheu IS away from the leads, but 
i'a.nnnins\r=R+w/2 = ''Piea.d whcu 9 IS in the range of a lead. 
Also, the radial derivative i9V'annuius/i9'' is equated with 
the longitudinal derivative dipieaA/dx when they meet at 
the outer arc of the annulus. The difference between the 
straight transverse cuts of the leads and the outer arcs of 
the annulus is neglected. The wave functions in the leads 
and the annulus are hence matched. Expanding the wave 
functions in different regions with sufficient numbers of 
modeS)2^ one gets a set of equations relating the coef- 
ficients of the modes in different regions. With a given 
energy E and specified in-going subbands, one can obtain 
the transmission probabilities in the out-going subbands. 

In the QID case, we will consider that the particle is in- 
cident from one lead, and its energies is below the second 
subband and hence the particle propagates only within 
the first subband. The resulting transmission probabil- 
ity is to be compared with that in the ID case. We 
will use the more convenient longitudinal wave number 
k\\ = ^2m(E — i?Q'™-)/S instead of the energy i?, where 
i^o""- is the energy of the nodeless ground state of the 
isolated annulus in an individual case^ Defining a lon- 
gitudinal wavelength A|| by A|| = 27r/fc|| implies that the 
long- wavelength limit we consider is at An ^ w. Here we 
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FIG. 3: (Color online) The transmission probability T is 
plotted versus the dimensionless longitudinal wave numbers 
(2L/A|| in the case of QID channels, and 2L/A in the case 
of ID channels), for the case of Li = L2 = L. The QID 
results (solid lines), ID results due to our scheme [dashed 
iy = 1.9) and dash-dotted {v — 1) lines], and ID results 
due to the Grifflth scheme (dotted lines) are shown. T is 
plotted for the cases of broad and narrow QID channels, (a) 
R/w = 3.5 and (b) R/w = 9.5. Note that the Griffith result 
is independent of the channel widths. In the narrow channel 
case [(b)], the difference between the QID and u — 1.9 results 
is indiscernible in the scale of this graph. 



define the arm lengths by Li_2 = Rdi,2, where 61^2 are 
the angles shown in Fig.|2b). 



III. COMPARISON BETWEEN RESULTS 

Figure 121 shows the transmission probabilities obtained 
by different schemes, for the case of symmetrical arms 
in the ring {Li ~ L2 = L). We have considered broad 
[Fig. OJa)] and narrow [Fig. Ofb)] channels, and in both 
cases we have presented the result of the QID calculation, 
the result of the Griffith scheme, and the results of our 
scheme (with v = 1 and 1.9). The Griffith result is seen 
to differ very much from the QID result in all cases. The 
1/ = 1 scheme qualitatively captures the trend of change 
in the QID result when the channel width is changed, 
while the v ~ 1.9 scheme captures the QID result most 
satisfactorily. 

Besides the Griffith result, it is seen that all results 
in Fig. 131 show Breit-Wigner (BW) resonance peaks 
These BW peaks become sharper and shift toward the 
left when the channels narrow down, i.e., w/R [com- 
pare Figs. |3{a) and (b)]. Those peaks are due to the 
quasibound levels in the arms, and they are seen to be 



always blue shifted2» from the exact levels. In the ID 
case, the exact levels are at 2L/A = kL/ir = integer. 
The quasibound levels and the blue shift are results of 
the presence of an attractive potential at a Y-junction^ 
While the attractive potential in our scheme is mani- 
fest [see Eq. (0], the potential at a Y-junction of QID 
channels is not so obvious, but can have an intuitive un- 
derstanding as follows. Since a particle feels less con- 
fined at near a junction, the "band-bottom" at the vicin- 
ity of a junction is effectively lower, and therefore the 
region acts as an attraction center. This potential be- 
comes stronger when the channels become narrower, and 
that leads to the sharper and less blue shifted BW peaks 
[compare Figs.|3Ia) and (b)]. The growth of the potential 
at narrowing channels can be understood as a result of 
the departure from the case of very broad channels (i.e., 
w ^ R), in which the system has no difference in the 
"band-bottom" everywhere. 

While the result from the Griffith scheme is indepen- 
dent of the channel width and disagrees with the QID 
result, our scheme captures the trend of change in the 
transmission probability when the channel width is var- 
ied. Therefore, our scheme has correctly included the 
attractive nature of the clean Y-junction of QID chan- 
nels, though the strength has been underestimated (i.e., 
v ~ 1.9 is prefered to 1/ = 1). The misjudgment of 
an appropriate value for the parameter v is due to the 
fact that the details of the shape of the Y-junctions of 
QID channels and the actual dimensionality of the chan- 
nels are relevant. For instance, our simple TB argument 
which leads to Eqs. ((Sj and lO does not show the differ- 
ence between junctions with different relative directions 
of branching channels, and also does not distinguish a 
three-dimensional (3D) cylinder from a 2D strip as a QID 
channel. But in reality, the appropriate parameter i/'s in 
those different cases may likely be different. In the 2D 
cases we have just seen in Fig. O the same kind of Y- 
junction has been involved, and the effective potential in 
our scheme is characterized by an almost constant 1/ in 
both the broad [Fig.l^Ja)] and narrow [Fig.|2|[b)] channel 
cases. 

Therefore, though the parameter 1/ can not be derived 
analytically, it can be readily fixed for a particular kind of 
junction by comparing the ID result with the QID result. 
What we have done in Fig.|31has been a comparison which 
involves a tedious calculation. Actually, other simpler 
comparisons also work. For instance, one may consider 
the bound state at the junction due to the attractioni^ 
On one hand, for a junction of three ID channels, with 
the channels extended to infinity like what we depict in 
Fig.^c), the negatively valued bound state energy E can 
be readily found by using ipi = e~'^^\ where i = 1, 2, 3, 
K = \/—2mE/h, and Eq. The energy E is found to 
be lower than zero by an amount of 2v^h^ / {Qmw^). On 
the other hand, the bound state at a T-shaped junction 
of three QID channels, with the channels extended to 
infinity like what we depict in Fig.^a), was studied by 
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FIG. 4: (Color online) The transmission probability T is 
plotted versus the dimensionless longitudinal wave numbers 
[(Li + L2)/A|| in the case of QID channels, and (Li + I/2)/A 
in the case of ID channels], for the case of asymmetrical arm 
lengths. The QID results (solid lines), ID results due to our 
scheme [dashed {v = 1.9) and dash-dotted {u — 1) lines], 
and ID results due to the Griffith scheme (dotted lines) are 
shown. T is plotted for the cases of broad and narrow QID 
channels, with small and appreciable differences in the arm 
lengths. T is plotted for (a) R/w = 3.5, L2/L1 = 0.9, (b) 
R/w = 9.5, L2/L1 = 0.9, (c) R/w = 3.5, L2/L1 = 0.7, 
and (d) R/w = 9.5, L2/L1 = 0.7. Note that the Griffith 
result is independent of the channel widths. In the narrow 
channel cases [(b) and (d)], the differences between the QID 
and 1/ = 1.9 results are indiscernible in the scale of this graph. 



Schult et alm^ The energy of the state was numerically 
found to be lower than the first subband bottom by an 
amount of QAQn^h^ /{2mw^)^ Equating the two ener- 
gies in the ID and QID cases, one gets i' ~ 2.05, which 
is about the number we use in Fig. O and as we will sec, 
that in Fig. H 

Figure 21 shows the transmission probabilities for the 
case of asymmetrical arm lengths. It is seen that in 
all cases, broad channels [Figs. El^a) and (c)] and nar- 



row channels [Figs.^fb) and (d)], small difference in arm 
lengths [Figs.^fa) and (b)] and appreciable difference in 
arm lengths [Figs.^c) and (d)], there are good compar- 
isons between the results due to our h' = 1.9 scheme and 
the QID calculation. All the essential features, such as 
the relative positions of the BW and Fano profiles^S in 
the QID results, are nicely reproduced. Note that the 
number 1.9 agrees with the one used in Fig. |31 

For the ID models, including Griffith's and ours, the 
perfectly zero transmission dips of the Fano profiles 
are located exactly at the cigenenergies of an isolated 
ringffi with a circumference of Li + L2. In the ID 
case, these cigenenergies are exactly at (Li -I- L2)/\ = 
k{Li + L2)/{2tt) = integer. In the QID cases shown in 
Figs. 13 and ^ the cigenenergies are numerically found 
to be at the [Li -I- L2)/A|['s deviated by not more than 
0.5% from the integers on the horizontal axes. For the 
QID model, we find that those zero transmission dips 
may coincide with the cigenenergies of an isolated annu- 
lus only in the long- wavelength limit. As in the case of 
symmetrical arms, the Griffith result disagrees with the 
QID result, and our simple TB argument which leads to 
Eqs. lO and © has underestimated the strength of the 
effective potential at the junction, i.e., v = 1.9 is prefered 
to = 1. 



IV. CONCLUDING REMARKS 

It is seen that in all the above cases the results from 
the Griffith scheme are not in congruence with the QID 
results. The Griffith result is regardless of the channel 
width, whereas the QID result shows a strong depen- 
dence on that. Our model gives a result in much better 
agreement with the QID result. The trend of change in 
the transmission probability and the relative positions of 
the resonance profiles are impressively reproduced. In 
view of these calculations, it is clear that the Griffith 
scheme which is frequently adopted in the literature, does 
not describe a clean junction of QID channels and is 
definitely not for the ID limit of the QID models. In 
the small width limit, a Y-junction of QID channels is a 
strong scatterer, and that makes the QID system stud- 
ied in this paper not at all an "open" system. Speaking 
reversely, adding a repulsive potential to a Y-junction of 
QID channels may weaken the scattering effect and en- 
hance the transmission through the junction at low ener- 
gies, and away from the levels. When a strong magnetic 
field is present, our model may not apply since the field 
creates an additional asymmetric transverse confinement. 

In conclusion, we have proposed a connection scheme 
with a parameter at a Y-junction of ID channels. The 
parameter v can be most easily fixed by comparing the 
energy of the bound state at a Y-junction of QID chan- 
nels, to the energy of the bound state at a Y-junction 
of ID channels due to Eq. ©. The scheme refiects the 
presence of an effectively attractive potential at a clean 
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Y-junction of QID channels. The disregard of this po- 
tential in the Griffith scheme makes its result compares 
poorly with the exact QID result. 
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